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Abstract The U(5) — 0(6) transitional behavior of the Interacting Boson Model in the large N 
limit is revisited. Some low-lying energy levels, overlaps of the ground state wavefunctions, B(E2) 
transition rate for the decay of the first excited energy level to the ground state, and the order 
parameters are ealeulated for different total numbers of bosons. The results show that eritieal 
behaviors of these quantities are greatly enhanced with inereasing of the total number of bosons N, 
especially fractional occupation probability for d bosons in the ground state, the difference between 
the expectation value of in the first excited O"'' state and the ground state, and another quantity 
related to the isomer shift behave similarly in both the 0(6) — U (5) large N and U (5) — SU (3) phase 
transitions. 
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Quantum phase transitions have been attracting a lot of attention in many areas of physics. This 
is understandable because they are very important for gaining a deeper understanding of various 
quantum many-body systems. [^“^1 In atomic nuclei, such quantum phase transitions can be related 
to different geometrical shapes of the system, which can be described either by the Bohr-Mottelson 
modelW (BMM) or by the Interacting Boson Model^^l (IBM). As summarized by lachello,!®! the 
study of shape phase transitions in atomic nuclei was initiated in the early following some 

previous workl^l by Gilmore. It is now widely accepted that the three limiting cases of the IBM 
correspond to three different geometric shapes of nuclei, referred to as spherical (vibrational limit 
with U(5) symmetry), axially deformed (rotational limit with SU(3) symmetry), and y-soft (triaxial 
with 0(6) symmetry). This picture is captured by the so-called Casten triangle.More interesting 
scenarios occur when a system is in between two different phases, in which case a quantum phase 
transition occurs at the corresponding critical point. A critical point at finite N with E(5) symmetry 
along the 0(5) — 0(6) leg of the Casten triangle was shown to exist in [11], and many examples 
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confirming the nature of this transition in realistic nuclear system have been reportedJ^^l Recently 
this transitional region has been studied for relatively large N values and the results show that the 
critical point region becomes progressively narrower as the boson number N increases.This 
phenomenon has been explained in [14] in terms of a quasidynamical symmetry. In order to study 
the large N limit situation corresponding to the classical BMM, one must approach the large N 
limit from results for finite N if algebraic results for the large N limit is not available. In this 
Letter we revisit the U{5) — 0(6) transitional case in the large N in detail to see whether there are 
substantial changes that occur as N grows ever larger, which serves as a supplement to the results 
reported in [13] and [14]. 

Our investigation is based on the following schematic 0(5) — 0(6) Hamiltonian: 

H = {l-x)h, + j^^S+S-, ( 1 ) 

where hd = J2 dl^dm is the total number of d-bosons, S~^ = ^(d^ • d^ — s^^) and S~ = ^(d • d — s^) 

m 

are generalized boson pair creation and annihilation operators, f{N) is a linear function of total 
number of bosons N, and x is the control parameter of the model. It should be obvious that the 
system is in the 0(5) limit when x = 0 and in the 0(6) limit when x = 1. As the control parameter 
X varies continuously within the closed interval [0, 1], the system described by (1) undergoes a 
shape (phase) transition from 0(5) to 0(6). 

To diagonalize Hamiltonian (1), we expand the eigenstates of (1) in terms of the 0(6) D 0(5) D 
0(5) D 0(3) basis vectors \N Ud v LM) as 

\N e u OM; x) = ^ 0|^(x)|iV Ud v LM), (2) 

where 0|^(x) is the expansion coefficient, ^ is an additional quantum number needed to label 
different eigenstates with the same quantum numbers v, L, and, M. 

To show how the energy levels change as a function of the control parameter x and the total 
number of bosons N, the lowest 25 energy levels as a function of x for a system with fixed quantum 
number u and f{N) = N for N = 10, 40, 120, and 300 are shown in Fig. 1. It can be seen from 
these results that there is a minimum in the low-lying excitation energy when the control parameter 
has a value in the range 0.45 < x < 0.65, with the minimum growing sharper as the total number 
of bosons increases. This control parameter region is recognized as the critical point region of the 
vibrational-gamma soft transition. To the left of the critical point region, 0 < x < 0.45, there are 
9 degenerate levels (x = 0) that gradually split with increasing x into 25 non-degenerate levels. 
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Fig 1. The lowest 25 energy levels (in arbitrary unit) of Hamiltonian (1) with f{N) = N as & function of x for N = 10, 40, 120, 


and 300, respectively. 






Fig. 2. Overlaps of the ground state wavefunction, where the full line shows the overlap 1(0^; ir|0g; a: = 0)|, and the dotted line 
shows the overlap \{^g\x\xg\x = 1)|. 

Similarly, beyond the critical region, 0.65 < re < 1, the 25 non-degenerate levels coalesce into 5 
degenerate levels (x = 1). It should be noted that apart from the end points, the levels are truly 
non-degenerate, and that the level density grows rather dramatically within the critical point region 
with increasing N] and furthermore, as N grows the critical point region becomes progressively 
narrower with a cusp around x ~ 0.45, which is in agreement with the observation reported in [14], 
in which only the V = 40 case was shown. 

The corresponding overlaps of the ground state wavefunctions of Hamiltonian (1) as a function 
of the control parameter x with those of limiting cases 1(0^; xjO^; xq) | for xq = 0, 1 for different 
total number of bosons N = 10, 40, 300, 1000 were also calculated and the results are shown in 
Fig. 2. It can be seen from Fig. 2 that there is a cross-over point at a certain nonzero amplitude 
around x ~ 0.65 for the overlaps 1(0^; x|0g; x = 0) | and 1(0^; xjx^; x = 1) | when N is relatively small, 
which yields to a cross-over region with near zero amplitude when N becomes larger. Furthermore, 
there is a sharp change in |(0g;x|0g;x = 0)| around a critical point Xc ~ 0.45 in the large N limit. 
These results suggest that the largest absolute value of the derivative of |(0g;x|0g;x = 0)| with 
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Fig. 3. B{E2) transition rates for decay of the first excited 1 ^ = 1 energy level to the ground state for N = 10, 40, 300, and 
1000 expressed in units with B{E2; 1 —> 0) = 100 in the f/(5) (a^ = 0) limit. 

respect to x occurs around the critical point in the large N limit. While both |(0g;x|0g;x = 0)| 
and |(0g;x|0g;x = 1)| are all rather smooth in the relatively small N cases. 

B{E2) transition rates for decay of the first excited u = 1 energy level to the ground state for 
N = 10, 40, 300, and 1000 expressed in units with B{E2; 1^0) = 100 in the 17(5) {x = 0) limit 
were also calculated. The E2 transition operator was chosen as T{E2) = e 2 {s^d + d)s)'^, where 
62 is the effective charge. The results are shown in Fig. 3. It can be seen quite clearly that the 
B{E2; 1^0) changes rather smoothly with x for small N, while there is a sharp change at the 
critical point when N is large enough. This behavior of the B{E2; 1^0) was also reported for the 
V < 60 cases considered in [14]. 

The fractional occupation probability for d bosons in the ground state, pd = {nd)/N as a 
function of x was reported in [13] and [15]. It was shown that an order parameter to signify a 
second-order phase transition can be chosen to be pd- Our calculation indicates that the system 
is almost in the 17(5) limit when x ~ 0 — 0.45 in the large N limit, which corresponds to an 
s-boson condensate. The occupation probability pd gradually increases within the critical region 
for relatively small N values with the change in pd becoming sharper and sharper with increasing 
N, which is in agreement to the results reported in [13] and [15]. Since the behavior of the order 
parameter pd is the same for both first- and second-order transitions for the small N cases, in order 
to distinguish whether the phase transition is of first or second order from model calculations, 
another order parameter, the difference between the expectation value of in the first excited 0 "*“ 
state and the ground state, vi = Q!o(( 02 |nrf| 02 ) — (0g|n(i|0g)), was introduced in [15]. The authors 
showed that vi displays a wiggling, sign-change-behavior in the region of the critical point due to 
the switching of the two coexisting phases, which is characteristic of a first-order transition, while 
vi has a smoother behavior that is characteristic of a second-order transition. It should be pointed 
out that the conclusion made in [15] are for finite N only. However, since order of a phase transition 
should always be defined in the thermodynamic limit, an effective order parameter must also behave 
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Fig. 4. Behavior of the order parameters vi and V 2 as functions of the control parameter x for different N values, where the 
parameters oq 00 and V 2 , respectively, are set to be 1. 


differently in phase transitions with different orders. To see whether the order parameter vi and 
another quantity V 2 = /3o((2i|nd|2i) — (0i|n(i|0i)) related to the isomer shift S{r‘^) = {r‘^) 2 i — (^^)og 
introduced in [15] satisfy this criterion, both vi and V 2 were calculated for the N = 10, 40, 300, 
and 1000 cases. The results are shown in Fig. 4. In order to compare curves of vi and V 2 for the 
different N cases, the parameters ao and Po were taken to be 1 instead of the 1/N used in [15]. 
Our calculation shows that; (a) both vi and V 2 have a smooth behavior when N is relatively small; 
(b) vi gradually displays of a sign-changing nature in the critical region when N is large enough, 
with this behavior being greatly enhanced in the large N limit; and (c) there is an obvious peak 
in V 2 in the large N limit, while V 2 is rather smooth for relatively small N. These results shown 
in Fig. 4, together with those shown in [15], indicate that the order parameters vi and V 2 , like 
another order parameter pd, behave similarly in both the 0(6) — U{5) large N and U{5) — SU{3) 
phase transitions. Due to current computation limitation, one can not calculate these quantities 
in the U (5) — SU (3) transitional case exactly for N > 30. Therefore, whether these quantities in 
the U (5) — SU (3) case will change substantially in the large N limit is still an open question. For 
relatively small N cases, however, as indicated in [15], the order parameters vi and V 2 are indeed 
behave differently in the 0(6) — 0(5) and 0(5) — 5*0(3) phase transitions, which, therefore, can be 
used to signify the order of the transition from the small N cases. 

In summary, the 0(6) — 0(5) transitional behavior in the large N limit has been revisited. 
Some low-lying energy levels, overlaps of the ground state wavefunctions, B{E2) transition rates 
for decay of the first excited u = \ energy level to the ground state, and the order parameters vi 
and V 2 related to the isomer shifts were calculated for different total number of bosons. It is found 
that the critical behaviors of these quantities are greatly enhanced with increasing of the total 
number of bosons V, especially all the order parameters, pd, vi, and V 2 behave similarly in both 
the 0(6) — U{5) large N and U{5) — 50(3) phase transitions. The drastic enhancement of these 
quantities near the critical point may be explained in terms of a quasidynamical symmetry.The 
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“specific heat” introduced in [16] seems also suitable to be used to classify the order of the phase 
transitions since these quantities behave quite differently in first and second order phase transitions 
even when the N is finite. 

In the IBM for atomic nuclei the total number of bosons N is phenomelogically related to be 
the number of valence s and d nucleon pairs, which is usually a relatively small number. However, 
in the large N limit, the IBM yields to the BMM, in which there is no restriction on the number 
of bosons; indeed, in principle, this should correspond to the N ^ oo limit. Therefore, the results 
shown in this Letter should be helpful in understanding the nature of the vibration to gamma- 
soft phase transition in the BMM. It is interesting to check to see whether there are substantial 
differences between the E{5) symmetry derived from an extreme case of the BMM and systems 
described by a 17(5) —0(6) Hamiltonian with the finite N based on the IBM. A recent study suggest 
that the E{5) symmetry can only be described approximately in the IBM,[^^] which is a conclusion 
that is consistent with our results. 
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